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ABSTRACT 



Six intensity groups with ~ 150 BATSE gamma-ray bursts each are 
compared using average emissivity curves. Time-stretch factors for each of the 
dimmer groups are estimated with respect to the brightest group, which serves 
as the reference, taking into account the systematics of counts-produced noise 
effects and choice statistics. A stretching/intensity anti-correlation is found with 
good statistical significance during the average back slopes of bursts. A stretch 
factor ~ 2 is found between the 150 dimmest bursts, with peak fiux < 0.45 ph 
cm^^ s^^, and the 147 brightest bursts, with peak fiux > 4.1 ph cm~^ s^^. On 
the other hand, while a trend of increasing stretching factor may exist for rise 
fronts for burst with decreasing peak fiux from > 4.1 ph cm~^ s^^ down to 0.7 
ph cm^^ s^^, the magnitude of the stretching factor is less than ~ 1.4 and is 
therefore inconsistent with stretching factor of back slope. 
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1. Introduction 

Gamma-ray bursts are known to have very different intensities. The peak ffux of 
bursts varies more than 2 orders of magnitude from the BATSE trigger threshold of about 
0.3 ph cm~^ s~^ up to the highest measured values of about 50 ph cm"^ s^^ ( [Fishman] 
and Meegan 1995|) . The rich statistics of the BATSE 4B Catalog ( Paciesas et al. 1998| ) 



enables one to divide all detected bursts into several intensity groups with reasonably large 
numbers of events (~ 150) in each of them. Although individual bursts have very different 
time profiles, a characteristic average temporal signature can be produced for each group. 
Using these average signatures, a comparison can be made between the duration of different 
brightness groups of gamma-ray bursts. The main goal of this comparison is to test the 
brightness-dependent stretching of gamma-ray bursts. 

It is well-known that individual pulses from radio pulsars have quite variable emission 
time profiles, so that the characteristic periodic pulsar signal is hardly recognizable by 
examining a short interval of real-time data. The randomized time profile is transformed 
into the stable signature of a pulsar light curve by using the epoch folding technique, 
averaging the data over many periods. A similar signature has been proposed for cosmic 
gamma-ray bursts: averaging the time profiles of individual events by the normalized 
peak-alignment technique, where each time profile is normalized by the peak number of 
counts Cjnax; aligned at peak time bins tmax and then averaged for all bins along the time 
scale (Mitrofanov et al. 1996, hereafter [Paper 1] ). This technique produces an Average 
Curve of Emissivity, or ACE, and the corresponding signature represents the averaging of 
observed time histories of bursts. 

ACE curves have already been studied for different energy ranges and for different 
intensity groups of bursts, and they have been found to be rather convenient signatures to 
describe the basic properties of the slow temporal variation of bursts ([Mitrofanov et ar 
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1997). This signature averages out the fast variations of bursts in the sample and describes 
the general envelope of rising emission before the main peaks of bursts and the subsequent 
decaying tail. An analytic approximation to the ACE has been found, which results in a 
very acceptable fit to the observations. The ACE equivalent width parameter describes the 
mean time scale of the slow variation of GRBs. 

The comparison of the generic time-dependent properties of gamma-ray bursts with 
different brightnesses has recently become an issue of common interest. Indeed, if the 
difference in intensities represents the difference in distances to the emitters, then the 
difference between the average time signatures of different brightness groups should 
manifest effects of cosmological time-dilation, where dimmer bursts are observed to be 
broader than brighter events because they were emitted at larger cosmological distances. 
The time-dilation test based on the comparison of the averaged temporal structures for 
dim and for bright sets has already been done by two groups, which have drawn two 
opposite conclusions: in the first case a large time-dilation effect ~ 2 was seen ( [Norris et al 



19941 , [Norris et al. 199"7D , in the other case none was observed ( [Mitrofanov et al. 1994| & 



Paper H ). 



The effect of cosmological time-dilation should also be accompanied by the effect of 
photon energy red-shift, which could influence the time-stretching value ( [Fenimore fc Bloom 



|1995a| , [Mitrofanov et al. 1997|) . On the other hand, sources of bright and dim gamma-ray 



bursts could be non-standard candles with intrinsic luminosity-based correlations ([Brainerd 



[1997[) . In this case, the observed time-stretching of dim bursts could result from the fact 
that shorter outbursts correlate with smaller intrinsic luminosity of their sources. Also, 
one should take into account possible distance-related evolution of bursts emitters: in the 
co-moving reference frames, more distant sources could have shorter time profiles than 
nearer sources. These intrinsic effects could interfere with the cosmological effects. 
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Therefore, in order to make conclusive statements about generic time-dependent 
properties of gamma-ray bursts, one should distinguish between the physical effect of time- 
dilation, which results from the expansion of the Universe, and the phenomenological effect 
of time-stretching, which is found by the comparison of bursts with different intensities. 
We think that the first logical step should be the investigation of the stretching/intensity 
phenomenon by comparing time profiles for different intensity groups of bursts. If a 
significant time-stretching effect is found between different intensity groups, then the 
next logical step can be made. This step should be the physical interpretation of the 
observed stretching, according to predictions of a cosmological model, taking into account 
the geometrical effects of time-dilation and red-shift, the physical effects of time-energy 
dependence of emission, the effects of a broad luminosity distribution of emitters, and 
finally, the effects of distance-related evolution. In this paper, we restrict ourself to purely 
phenomenological studies of the stretching/intensity correlation of bursts, using the average 
emissivity curves for different brightness groups of bursts as the tool for measuring the 
stretching. 



A total set of 887 bursts with tgo > 2 s from the BATSE 4B Catalog (|Meegan et al 



1994, Faciesas et al. l^UE) have data available on the 1024 ms timescale. We have divided 



these into 6 intensity groups (Table 1), using the peak flux parameter F^^^^^ as the selection 
criterion. For each of these groups, an ACE proflle is produced using counts observed in the 
broad energy range 50-300 keV. The technique used to produce ACE proflles is described 
in 



Paper I| . Below, the stretching/intensity effect is investigated and stretching factors are 
evaluated between each ACEj for the dimmer groups (i = 2 - 6) and the reference group 
(ACEi). We have done this analysis for the total proflle as well as for the rise fronts and 
the back slopes separately. 
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2. Analytical approximation of ACE profiles 

We have previously determined that there is a simple analytical form that fits the 
ACE profiles of all burst intensity groups in several different spectral ranges quite well 
( [Mitrofanov et al. 1997|) . For each selected group (i) of bursts, the analytical approximation 



to the ACE is 






with different power indices (a^p, agg) at the rise front (RF, with t < tmax) and at the back 
slope (BS, with t > tmax), respectively. 

Another approximation for ACE-like profiles with exponential wings has been suggested 
by [Stern et al. 1996| . A direct comparison between these two approximations has shown 
that a power law (eq. 1) fits the 1024 ms time scale ACE profiles over the range of 20 time 
bins before and after the peak much better than the exponential model. For example, with 
intensity group 1 (Table 1) we found a reduced x^ = 1-51 for an exponential-type law, while 
Eq. 1 provided a better fit, with reduced x^ = 0.88. For this reason, we use Eq. 1 to model 
ACE profiles in this paper. 

3. Two procedures to estimate stretch factors between ACE curves of 

different intensity groups 

When cosmological models were first suggested for GRBs, the observational tests for 
cosmological signatures were expected be quite obvious; for the case of ACE stretching 
coefficients, a factor of at least 2 should be observed between the brightest and dimmest 
burst intensity groups. Therefore, a simple procedure was devised to look for evidence for 
this stretching (see [Paper i| ). However, a stretching factor of the expected magnitude was 
not found; indeed, recently-developed cosmological models for bursts predict much smaller 
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average stretching, despite the fact that some emitters should exist at large cosmological 
red-shifts (see e.g. [Brainerd 1997|) . In order to make further progress, a more sensitive 



comparison method between different ACE profiles has been developed, allowing the 
measurement of weak-signature stretching factors (less than 1.5). To be successful, this 
method must take into account several statistical biases associated with the limited number 
of bursts in different intensity group samples and with the poor counts statistics of weak 
events. 

We have developed two such procedures to make a robust comparison of ACE 
profiles between different brightness groups of bursts. Procedure (a), which is relatively 
time-consuming, requires the following steps to compare intensity groups i = 2 - 6 with the 
reference group i = 1: 

(al) An arbitrary stretching coefficient Y is selected, and each burst of the reference 
group 1 is stretched by the factor Y. 

(a2) An artificially-dimmed reference group, which has the same signal-to-noise ratio 
as group i, is created from bursts of the stretched reference group (see Section 5). 

(a3) ACE profiles are produced for the stretched and dimmed version of the reference 
group and for the original group i. The value of x^ is used to measure the difference between 
them. To find the best fit stretching coefficient Yi, the value of the stretching coefficient Y 
is changed, and the cycle (al)-(a3) is performed again to minimize x^- 

The time-efficient procedure (b) is more simple and straight-forward than the first one. 
It consists of the following steps: 

(bl) An artificially-dimmed reference group that has equal signal-to-noise ratio to the 
the group i is created from the reference group (see Section 5). 

(b2) An ACE profile is produced for the artificially-dimmed reference group, and the 



parameters of the best-fit analytical model (Eq. 1) are computed. 

(b3) The stretching parameter Y for the time constant At'-^-' is introduced into the 
analytical model (Eq. 1) for the artificially dimmed reference group 1: 

Fixing all other parameters, the best fit value for Y, as a free parameter, is found for the 
ACE of each group i = 2 — 6. This value is taken to be the stretch coefficient Yi between 
the group i and the reference group. 

A direct comparison of these two procedures has been performed for the back slopes 
(time bins after the peak in the ACE) of bursts using the 3B Catalog data base ( |Litvak| 
let al. 1997a ). Excellent agreement between the best-fit stretch factors (Y) estimated from 



both procedures was found. The difference between them is much less than la for choice 
statistics (see Section 4). Therefore, the computationally time-efficient procedure, (b), is 
used below in studies of ACE stretching. 

4. Choice statistics and the related errors of the stretching coefficients 

The proposed method suggests that we may use the ACE as the signature representing 
the slow variability of a selected sample of gamma-ray bursts. Individual bursts have very 
different time profiles, and one must use a very large number of individual events to build 
up a representative ACE that will be the same for different samples of bursts with similar 
brightnesses. 

The errors associated with ACE curves may be estimated from the sample variance 
for the selected groups of bursts. ACE profiles for two independent groups are statistically 
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indistinguishable, provided the difference between them is within these errors. The direct 
comparison of ACE profiles for different groups has shown that they are much more different 
than might expected from the variance within each sample only. This means that these 
groups are not representative samples, and a random selection of bursts within these groups 
does not ensure well-weighted contributions from all kinds of possible profiles. 

To study the random choice statistics of bursts, a Monte Carlo test has been performed 
using the total set of 603 bursts from the 3B Catalog QMeegan et al. 1997| ) that have 



catalog values of peak flux and durations T90 > 2 s (|Litvak et al. 1997b| ). Different numbers 



of bursts, A^ = 100, 150, 200, 250, 300, were assumed. We randomly selected 10^ times, 
two testing groups of A^ events each from the total set of bursts, produced ACEs for 
these groups, and then determined the best-fit stretch factor ychoice between them. The 
distributions of Ychoice for two different values for A^ are presented in Fig. 1. 

The spread of stretching coefficients ^choice due to random choice statistics is found 
to be much broader than would be expected from the sample variance for each individual 
group. Of course, the distribution of l^choice becomes narrower for larger A^, but even for the 
largest value, A^ = 300, it is still quite broad. 

Therefore, we conclude that for any two intensity groups, with A^ bursts each, the 
estimated stretching coefficient must be compared with standard deviations from the 
random choice distribution, evaluated for A^ events. For A^ = 150, the level of la significance 
corresponds to the stretching coefficient ~1.10. Therefore, for y ~ 1, one should use the 
errors 6Y ~ 0.10 for the la errors of the stretching coefficients. In the arbitrary case where 
A^ > 100, one can use following for estimations of errors: 



Below, these errors were used to estimate the significance of stretching between selected 
brightness groups (Tables 2 and 3). 
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5. Count-noise produced effects on the ACE 

The procedure to create each ACE includes the selection of the highest peak of each 
burst, with counting rate Cmaxj and the normalization of the time profile by the Cmax value. 
Therefore, the ACE is sensitive to a bias where positive fluctuations of counts (C + 6C) 
dominate the selection of Cmax- When the normalization is performed, a profile is lowered 
(Cmax + 5C)/Cmax timcs. The effect should be larger for dimmer bursts, where the ratio of 
(Cmax + 5C')/Cmax IS larger. 

To evaluate the effect, the reference group 1 (Table 1) has been transformed by the 
procedure of Monte Carlo noisification into 5 artificial reference groups, which have peak 
fiux distributions similar to the corresponding distributions for the 5 observed dim intensity 
groups, i =2-6 (Table 1). The Monte Carlo transformation procedure includes the following 
steps: 

a) For each original burst of the reference group (1), some counterpart event is 
randomly selected inside the testing group i. 

b) The ratio of peak fiuxes / = -FmaxZ-^mL is estimated between the peak fiuxes of the 
original burst and its weaker counterpart. 

c) The time profile of the original bright burst is divided by the factor /, and counts for 
an artificially-dimmed version of the original burst are simulated using Poisson statistics: 

Dj = ^ + Cb„ (4) 

where Csj, Cb,j are the signal and background counts accumulated during the jth time bin 
for a burst from the reference group. 

Using equation (2), stretching coefficients Fnoise have been estimated between the ACEi 
for the original reference group and ACEf''' for the artificial reference groups {i = 2-6). 
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One can see the largest effect is tliat tlie ACEg''* for tlie dimmest group {i = 6) must 
be broadened by a factor of 1.22 to equal the width of the ACEi for the original reference 
group. Therefore, if the observed stretching factor between ACEs for the actual observed 
dim group and the reference group is equal to 1, it should be interpreted as evidence for real 
stretching on the order of l^noise between the dim and reference groups, because the ACE for 
the dim group is known to be narrowed systematically due to larger relative Poisson noise. 
This effect could be ignored when testing for much larger stretching coefficients above 2.0 
( [Paper 1 ), but it should certainly be taken into account when the expected stretching might 



be as small as the noise-produced effect (see below). Both procedures (a) and (b) (Section 
3) take this effect into account. 

6. Stretching factors between the ACE for different intensity groups 

6.1. Stretching with respect to the reference group 

Using procedure (b) to compare between ACEs, stretching coefficients were estimated 
for each intensity group (Figure 2 and Table 2). The group of brightest events (1) is used 
as a reference and the corresponding stretching coefficient for this group was defined to be 
1. For each group, the stretching factors Y^J , Igg and Itot were estimated with respect 
to group 1 for the RF and BS wings and for the total ACE profiles, respectively (Figure 2 
and Table 2) . 

The rise front portions of bursts do not manifest any significant increase of stretching 
with decreasing intensity of bursts. The stretch factor for the rise front of group 5 with 
respect to the reference group is the only one that is above the 3a level (in this case, 
3cr = 1.34). On the other hand, at the back slopes the stretching factors leg are larger 
than 3(T already between groups 2 and 1, and increase up to ~2 with decreasing brightness 
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(Table 2). The noise-produced systematic narrowing of the ACE is taken into account for 
these values, because the test groups (i = 2 - 6) are compared with the correspondingly 
noisified reference groups. In the case of the comparison with the original reference group 
(1), the stretching would be less. The lack of stretching during the rising portion of the 
ACE of the dimmest group (6) may be an instrumental effect, resulting from the deficit of 
slow-rising events near the trigger threshold (e. g. [Higdon fc Lingenfelter, 1996| ). We know 



from the non-triggered burst sample of [Kommers et al. (1997)] that the threshold effect may 



change the estimate for the rise-front stretching factor only for the dimmest group 6 in our 
sample. 

When the total ACE profiles are compared, the stretching factors have intermediate 
values between the corresponding factors for rise fronts and back slopes. One can see that 
the back slope stretching factors for dimmer groups {i = 3-5) are ~ 2a larger than those 
associated with the rise front (Table 2). That is why one should study the stretching 
phenomena of bursts separately for the rise front and back slopes. 

6.2. Stretching with respect to the second brightest group 

With decreasing burst intensities, the largest increase of stretching between successive 
intensity groups happens between the brightest and the second brightest groups (Table 2 
and Figure 1). Of course, this could result from a random choice fluctuation: the number 
of bursts in each group is still not large enough to exclude this possibility. Larger burst 
samples from subsequent BATSE catalogs should be used to check the jump of stretching 
for the two intensity groups bounded by the flux of ~ 4.1 ph cm~^ s~^. The jump may 
completely disappear for larger statistics; or, more interestingly, it could be confirmed as a 
real phenomenological effect. 
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However, one has to be sure that this jump-hke effect between the test group {i = 2) 
and the reference {i = 1) does not result from some systematic effect in the comparison 
procedure. Indeed, the reference group has been transformed by adding noise in order to 
be to be compared with dimmer groups, and some unknown systematics in this procedure 
could result in a jump of stretching between the reference group and dimmer groups. 

To check this possibility, the same procedure used for the estimations of stretching 
factors has been reproduced; however, this time assuming the second bright group (2) to be 
the reference (Table 3). The stretching factors Y^J and Ygg for group 2 have been defined 
as 1.0 in this case. To facilitate comparison between stretching coefficients based on group 
1 and group 2, the re-normalized factors 1.19 ■ Y^p and 1.35 ■ Ygg are also presented in 
Table 3. There is very good agreement between stretching factors Y^^^ (Table 2) and Y^"^^ 
(Table 3), based on the reference groups 1 and 2, respectively. 

Therefore, one can exclude the possibility that some systematic effect takes place when 
the reference group is compared with dimmer ones. One may conclude that in the 4B 
catalog data set a large stretching factor of about 1.35 really exists between the back slopes 
of bursts of intensity groups 1 and 2, separated by a flux ~ 4.1 ph cm~^ s~^. 

6.3. Consistency between stretching factors for intensity groups of the 2B and 

4B database 



A comparison between ACE profiles was done previously ( [Paper 1| ) using the database 



of 338 events from the BATSE 2B catalog ([Meegan et al. 1994] ), with durations Tgo > 1 



s. These were divided at the peak flux value ~ 1 ph cm^^ s^^ into two intensity groups 
with 143 bright bursts and 179 dim bursts. The estimated equivalent time widths tETW 
(see Paper I ) of ACE rise fronts, back slopes and total profiles for these old samples are 
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presented in Table 4. The values labeled "Old samples" are taken from Paper 1| but the 



errors are estimated according to random choice statistics as described above (Section 4). 
No significant stretching is seen between the ACE for these groups. 

The 4B catalog allows us to do a similar analysis for new samples with as many as 480 
bright bursts and 464 dim events, and to check the consistency between the estimations 
of equivalent time widths for the old and new samples. One concludes from Table 4 that 
the "old" and "new" results are consistent for intensity groups divided by the same peak 
fiux value ~ 1 ph cm^^ s^^. There is no evidence for stretching for the "old" groups with 
143 bright and 179 dim events either during the rise front or the back slope. Using the 
much better statistics available with the "new" intensity groups, there still is no significant 
stretching during the rise fronts, but there is some effect when r^TW ^^^ compared. 

The separation of all bursts into bright and dim groups by the peak fiux ~ 1 ph cm~^ 
s~^ was appropriate when the goal was to test for an obvious signal of stretching by factors 
large as ~ 2. On the other hand, the 2B catalog sample was too small for more accurate 
sampling. Now, with much better statistics, more subtle stretching effects can be seen for 
several intensity groups (Table 2), in particular at the back slope of ACE profiles. In the 
half-to-half separation, a stretching effect of ~ 2a is also seen at the back slope (Table 4). 

7. Conclusions 

Six burst intensity groups with ~ 150 events each have been compared to determine 
a stretching effect between the ACE profiles for dim and bright groups. To study the 
stretching, a separate comparison is preferable for the average rise fronts and back slopes. 
The Pearson x^ statistic between the reference and the stretched dimmer groups allows 
us to find the most probable stretching factors for the ACE of different intensity groups 
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(Table 2), and the resulting reduced x^ values are significantly smaller when the rise fronts 
and back slopes are each compared separately. There is a significant difference between the 
corresponding stretching factors for the rise fronts and back slopes in different intensity 
groups. During the rise front a stretch factor of 1.39 ± 0.14 is found between intensity 
groups 5 and 1, which seems to be marginally significant. The other factors determined for 
the rise fronts, taken together, also indicate some stretching effect, but not significantly. 

This study is based on BATSE 1 s resolution discriminator data. Higher temporal 
resolution datatypes begin either at the trigger time or 2 s before. Our choice of the 
continuously available 1 s data avoids possible systematic effects in the rise front due to 
mixing datatypes of differing temporal resolutions. 

There may be an instrumental triggering bias that selects against those dim bursts that 
rise more slowly on the average. This has been corroborated in the study of non-triggered 
bursts by P:<.ommers et al. (1997)| , who found that only some of the dimmest bursts failed 



to trigger due to slow rise times. The intensities of these events correspond to our dimmest 
intensity group 6. Since only group 6 is incomplete due to missing slow risers, and since 
this incompleteness affects only the rise front results, we base our conclusions on rise front 
stretching on groups 1 to 5. Therefore, omitting the dimmest group, there is the indication 
of a trend for the rise-front stretching factors (Figure 3). 

At the back slopes the estimated stretch factors are quite significant for all dimmer 
groups (2 = 2-6) with respect to the brightest one {i = 1), and the largest factor between 
the ACE profiles of dimmest 150 bursts and the brightest 147 bursts is about 2.1± 0.2. The 
non- stretching hypothesis may be significantly rejected for the back slopes of bursts. Our 
results are in qualitative agreement with recent estimations of time-stretching by [N orris et 



al. 1994| , where a different energy range was used for burst averaging and a different time 



scale was applied for burst selection. 
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We conclude that models of gamma-ray bursts should explain two phenomenological 
results: 

1) Back slopes for the 150 dimmest bursts, with peak fluxes < 0.45 ph cm~^ s~^, are 
on the average ~2 times longer in duration than the 150 brightest bursts, with peak fluxes 
> 4.1 ph cm^^ s^^. 

2) While a trend of increasing stretching factor, relative to group 1, may exist for rise 
fronts for groups 2 to 5, the magnitude of the stretching factor is less than ~ 1.4 and is 
therefore inconsistent with the stretching factor of the back slope. 

3)Finally, there is definitely no full-profile stretching between dim and bright bursts as 
large as ~ 2 or more. Whatever stretching exists is weak and so the correct determination 
of the stretching factor will require careful treatment of statistical and systematic effects 
that are comparable with the physical effect of stretching that we are looking for. Using the 
4B data, we conclude that a significant stretching by a factor of ~2 between the different 
brightness groups may only be resolved for the back slopes of the average light curves of 
GRB time profiles. 

Our primary concern has been assessing the observational evidence for time stretching 
in the average profiles of GRBs. The observation of stretching in the back slopes but not 
in the rise fronts of GRBs cannot be solely caused by cosmological effects - at least one 
of these phenomena seems to be intrinsic to GRBs, perhaps indicating a slope correlation 
with absolute intensity or the existence of source evolution. Whatever the explanation, 
determining the distance scale of GRBs using cosmological tests is proving to be more 
difficult than had been hoped. 
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Fig. 1. — The ychoice distribution for N=150 events is shown by a thick hne, while the 
distribution for N=300 events is shown by a thin hne. 

Fig. 2. — The ACE for the reference group {thick line) is compared with the ACE for dim 
groups {thin lines). 

Fig. 3. — The best-fit stretching factors for the average back slopes {solid) and rise fronts 
{dashed) relative to the brightest group {solid arrow) as a function of peak flux. The value 
for the rise front for the dimmest group 6 may be altered by missing "slow riser" triggers - 
see text. 
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Table 1. 



Intensity 
group 



Peak flux 



(ph cm 



-2 c,-l^ 



Number 

of bursts 



Y ■ 



1 
2 
3 

4 
5 
6 



>4.10 
1.76-4.10 
1.05-1.76 
0.67-1.05 
0.45-0.67 

<0.45 



147 
153 
148 
148 
148 
150 



1.00 
1.08 
1.11 
1.15 
1.19 
1.22 



Table 2. 



Intensity group 



Y, 



(1) 



RF 



Relative Stretching coefficients 



K 



(1) 



BS 



Fn 



(1) 



TOT 



1 


1.00 


2 


1.19±0.12 


3 


1.23±0.12 


4 


1.22±0.12 


5 


1.39±0.14 


6 


1.07±0.11 



1.00 
1.35±0.14 
1.54±0.15 
1.48±0.15 
1.78±0.18 
2.10±0.20 



1.00 
1.28±0.13 

1.44±0.14 
1.39±0.14 
1.60±0.16 
1.50±0.15 
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Table 3. 



Intensity group 




ri? ■ 1.19 


-•^BS 


Y^J ■ 1.35 


2 


1.00 


1.19 


1.00 


1.35 


3 


1.06±0.11 


1.26 
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Table 4. 



Parameters to compare 
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